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^vq We report the discovery of spatially-extended, non-thermal-like emission components 

in Chandra X-ray spectra for five of a sample of seven massive, merging galaxy clus- 

^ tcrs with powerful radio halos. The emission components can be fitted by power-law 

models with mean photon indices in the range 1.5 < F < 2.0. A control sample of 
\^ regular, dynamically relaxed clusters, without radio halos but with comparable mean 

thermal temperatures and luminosities, shows no compelling evidence for similar com- 
ponents. Detailed X-ray spectral mapping reveals the complex thermodynamic states 
of the radio halo clusters. Our deepest observations, of the Bullet Cluster IE 0657-56, 
demonstrate a spatial correlation between the strongest power-law X-ray emission, 
highest thermal pressure, and brightest 1.34GHz radio halo emission in this cluster. 
O We confirm the presence of a shock front in the IE 0657-56 and report the discov- 

^ ery of a new, large-scale shock front in Abell 2219. We explore possible origins for the 

C/3 power-law X-ray components. These include inverse Compton scattering of cosmic mi- 

I crowave background photons by relativistic electrons in the clusters; brcmsstrahlung 

from supra-thermal electrons energized by Coulomb collisions with an energetic, non- 
thermal proton population; and synchrotron emission associated with ultra-relativistic 
t> electrons. Interestingly, we show that the power-law signatures may also be due to 

complex temperature and/or metallicity structure in clusters particularly in the pres- 
ence of metallicity gradients. In this case, an important distinguishing characteristic 
OP between the radio halo clusters and control sample of predominantly cool-core clus- 

^"-^ ters is the relatively low central X-ray surface brightness of the former. Our results 

1— H have implications for previous discussions of soft excess X-ray emission from clusters 

and highlight the importance of further deep X-ray and radio mapping, coupled with 
OO new hard X-ray, 7-ray and TeV observations, for improving our understanding of the 

non-thermal particle populations in these systems. 

Key words: X-rays: galaxies: clusters - galaxies: clusters: individual: IE 0657-56, 
^ A665, A2163, A2255, A2319, A2744, A2219, A576, A1795, A2204, A478, A2029 - 

}^ radiation mechanisms: non-thermal - intergalactic medium - magnetic fields 
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1 INTRODUCTION 

Combined X-ray, optical and radio studies of massive, X- 
ray luminous galaxy clusters demonstrate a clear correla- 
tion between the presence of spatially extended radio halo 
emission and strong, recent merger activity {eg. Buote 2001; 
Schuecker et al. 2001; Govoni et al. 2001; Giovannini & Fer- 
etti 2002). The observed radio halo emission is synchrotron 
radiation from relativistic electrons, thought to have been 
re-accelerated by turbulence and/or shocks resulting from 
merger activity {eg. Ensslin et al. 1998; Bykov et al. 2000; 
Brunctti et al. 2001; Ensslin & Briiggcn 2002; Pctrosian 
2003; Brunetti et al. 2004; Brunetti & Blasi 2005; Feretti 



2005; Cassano, Brunetti & Sctti 2006; Cassano et al. 2007; 
Brunetti & Lazarian 2007; Pfrommer 2008; Pfrommer, En- 
ssUn & Springel 2008; Wolfe & Melia 2008; Petrosian & 
Bykov 2008). The surface brightness of the radio halo emis- 
sion is typically low (of order /xJy per arcsec^) but can span 
regions of a Mpc or more in size {eg. Feretti 2005; Feretti 
et al. 2005). 

The same population of relativistic electrons respon- 
sible for the radio halo emission is also expected to pro- 
duce non-thermal X-ray emission via Inverse Compton (IC) 
scattering of the Cosmic Microwave Background (CMB). 
In principle, the combination of radio and IC X-ray flux 
measurements can then be used to determine, or even map. 



© 0000 RAS 



2 E. T. Million and S. W. Allen 



the magnetic field strength in merging clusters [eg. Rephaeli 
1979; Fusco-Femiano et al. 1999; Eckert et al. 2008). How- 
ever, other processes can also contribute to a non-thermal- 
like X-ray emission signature]^ including shocks (which can 
temporarily heat regions of massive, merging clusters to 
temperatures of tens of keV e.g. Markevitch & Vikhlinin 
2007 and references therein; see also Komatsu et al. 2001; 
Allen et al. 2002), a range of non-thermal bremsstrahlung 
processes [eg. Ensslin, Lieu & Biermann 1999; Sarazin & 
Kempner 2000; Petrosian 2001; Liang et al. 2002; Dogiel 
et al. 2007; Wolfe & Melia 2008) and synchrotron emis- 
sion by ultra-relativistic electrons and positrons (Timokhin, 
Aharonian & Neronov 2004; Inoue, Aharonian & Sugiyama 
2005). Distinguishing such emission components in the pres- 
ence of potentially complicated astrophysical and instru- 
mental backgrounds, and the luminous, diffuse, ambient 
thermal X-ray emission from clusters, is challenging. 

At hard X-ray wavelengths, a detection of non-thermal 
X-ray emission from the nearby Coma Cluster has been 
reported by Fusco-Femiano et al. (1999; see also Fusco- 
Femiano et al. 2004) using BeppoSAX Phoswich Detection 
System (PDS) observations. The Coma cluster contains a 
powerful, extended radio halo and, being one of the near- 
est. X-ray brightest clusters, is an excellent candidate for 
the detection of IC X-ray emission. However, complexities 
associated with modelling the PDS data have led Rossetti 
& Molendi 2004 (see also Rossetti & Molendi 2007) to ques- 
tion the significance of this detection (but see also Fusco- 
Femiano, Landi & Orlandini 2007). Rephaeli, Gruber & 
Blanco (1999) and Rephaeh & Gruber (2002) also report 
a detection of non-thermal X-ray emission from the Coma 
cluster, based on Rossi X-ray Timing Explorer (RXTE) 
data. However, the large angular beam of the RXTE instru- 
ments also makes it difficult to disentangle true non-thermal 
emission components from spatial variations in the spec- 
trum of the diffuse intracluster gas. More recent observa- 
tions of the Coma Cluster with INTEGRAL, Swift/BAT 
and Suzaku, argue that there is no significant non-thermal 
X-ray emission at high energies in the Coma cluster (see Re- 
naud et al. 2006; Ajello et al. 2009; Wik et al. 2009). Eckert 
et al. (2008) report a detection of hard X-ray emission from 
the Ophiuchus cluster using observations made with the 
INTEGRAL satellite. However, more recent Swift/BAT 
and Suzaku analyses of the cluster do not confirm this result 
(Ajello et al. 2009; Fujita et al. 2008). Clearly the situation 
with regard to hard X-ray emission from clusters remains 
uncertain and more and better data are required. 

Soft X-ray (energy E < 0.5 keV) and far ultraviolet ob- 
servations offer another avenue to search for non-thermal X- 
ray emission. However, measurements at these wavelengths 
require precise background modelling and detailed knowl- 
edge of the thermal emission spectrum and Galactic column 
density along the line of sight {eg. Fabian 1996). Although 
strong arguments have been made for the presence of ex- 
cess UV/soft X-ray emission in the Coma Cluster {eg. Lieu 
et al. 1996; Bowyer, Berghofer & Korpela 1999; Arabad- 
jis & Bregman 2000; Bonamente et al. 2001a; Bonamente 



By non-thermal-like X-ray emission, we mean emission over 
and above that expected from the diffuse, ambient cluster gas 
that can be approximated by a power-law. 



et al. 2003; Finoguenov et al. 2003; Bowyer et al. 2004), for 
other systems, claims of diffuse UV emission over and above 
that expected from the intracluster plasma, remain contro- 
versial (see eg. Bregman & Lloyd-Davies 2006; see also Dur- 
ret et al. 2008). 

New and forthcoming 7-ray and hard X-ray missions 
such as the Fermt Gamma-ray Space Telescope (Fermi) 
and the Nuclear Spectroscopic Telescope Array (NuSTAR) 
should shed much light on the non-thermal and quasi- 
thermal particle populations in galaxy clusters. However, 
at present, arguably the best waveband in which to search 
for non-thermal and quasi-thermal emission signatures in 
galaxy clusters remains the ~ 0.5 — 8.0 keV energy range 
covered by major X-ray observatories such as Chandra and 
XMM-Newton. In particular, the Chandra X-ray Observa- 
tory offers excellent spatial resolution, good spectral resolu- 
tion, is well calibrated and has a well understood instrumen- 
tal background. The Chandra observing band extends down 
to energies soft enough to allow Galactic column densities 
to be measured precisely, and is wide enough that 'non- 
thermal' X-ray emission (be it true power-law-like emis- 
sion or complexities in the thermal emission spectrum) can, 
in principle, be detected. In particular, detailed spatially- 
resolved X-ray spectroscopy with Chandra minimizes sys- 
tematic uncertainties associated with spatial variations in 
the thermal spectrum of the diffuse cluster gas (although, 
as we show here, does not eliminate them). Such consider- 
ations have motivated the present study, in which we have 
used Chandra to carry out detailed spatially-resolved spec- 
tral mapping of the most X-ray luminous merging, radio halo 
clusters known, with the goal of searching for quasi-thermal 
and/or non-thermal X-ray emission. 

The structure of this paper is as follows. In Section 
2, we describe the target selection, data reduction, region 
selection, and spectral analysis. Section 3 summarizes the 
evidence for complex non-thermal-like X-ray emission com- 
ponents in the radio-halo clusters. Section 4 reports results 
from detailed thermodynamic mapping of the clusters, with 
particular emphasis paid to the two deepest observations, 
of the Bullet Cluster, IE 0657-56, and the relaxed, control 
cluster Abell 2029. Section 5 discusses the possible origins 
of the observed complex X-ray components. Section 6 sum- 
marizes our conclusions. A ACDM cosmological model with 
Ho = 70 km s~^ Mpc~^, fi^ = 0.3, and Qjy — 0.7 is assumed 
throughout. 

2 X-RAY OBSERVATIONS AND ANALYSIS 
2.1 Sample Selection 

Our targets are drawn from the sample of X-ray luminous 
clusters with large, extended radio halos discussed by Fer- 
etti & Giovannini (2007). These clusters have high qual- 
ity X-ray data available on the Chandra archive. All of the 
clusters show clear evidence for recent, major merger activ- 
ity and none of the clusters exhibits a cooling core (Fer- 
etti & Giovannini 2007). In order to minimize systematic 
uncertainties associated with modelling the background in 
the Chandra observations, we have imposed a lower redshift 
limit z > 0.03; this excludes the Coma Cluster from our 
study. We also require at least 2 x 10* net counts, which ex- 
cludes CL 0016+1609. Our target list of radio halo clusters 
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Table 1. Summary of the Chandra observations. Columns hst the target name, redshift, observation ID, detector, observation mode, 
exposure after cleaning and observation date. 



Name 


z 


Obs. ID 


Detector 


Mode 


Exposure (ks) 


Observation Da 








RADIO HALO CLUSTERS 






Abell 2319 


0.056 


3231 




V r AliN 1 


14.4 


Mar. 14 2002 


Abell 2255 


0.0809 


894 




r AilN 1 


38.2 


Oct. 20 2000 


Abell 665 


0.182 


3586 


A r^TG T 


V r AliN 1 


15.9 


Dec. 28 2002 


Abell 2163 


0.203 


1653 


A r^TG T 


V r AliN 1 


62.4 


Jun. 16 2001 


Abell 2219 


0.228 


896 


ACIS-S 


liAUN 1 


41.3 


Mar. 31 2000 


IE 0657-56 (1) 


0.297 


5356 


ACIS-I 


VFAINT 


96.2 


Aug. 11 2004 


IE 0657-56 (2) 




3184 


ACIS-I 


VFAINT 


82.4 


Jul. 12 2002 


IE 0657-56 (3) 


" 


5361 


ACIS-I 


VFAINT 


82.1 


Aug. 17 2004 


IE 0657-56 (4) 




5357 


ACIS-I 


VFAINT 


78.3 


Aug. 14 2004 


IE 0657-56 (5) 




4984 


ACIS-I 


VFAINT 


75.9 


Aug. 19 2004 


IE 0657-56 (6) 




4986 


ACIS-I 


VFAINT 


40.7 


Aug. 25 2004 


IE 0657-56 (7) 




5358 


ACIS-I 


VFAINT 


31.7 


Aug. 15 2004 


IE 0657-56 (8) 




4985 


ACIS-I 


VFAINT 


24.0 


Aug. 23 2004 


IE 0657-56 (9) 




5355 


ACIS-I 


VFAINT 


22.3 


Aug. 10 2004 


Abell 2744 


0.308 


2212 


ACIS-S 


VFAINT 


20.8 


Sep. 3 2001 


Abell 2744 




7915 


ACIS-I 


VFAINT 


14.5 


Nov. 11 2006 


Abell 2744 




8477 


ACIS-I 


VFAINT 


42.3 


Jun. 10 2007 


Abell 2744 




8557 


ACIS-I 


VFAINT 


23.5 


Jun. 14 2007 



NON-RADIO HALO CLUSTERS 



Abell 576 


0.0381 


3289 


ACIS-S 


VFAINT 


24.7 


Oct. 27 2002 


Abell 1795 (1) 


0.0634 


3666 


ACIS-S 


VFAINT 


13.2 


Jan. 18 2004 


Abell 1795 (2) 




5287 


ACIS-S 


VFAINT 


14.3 


Jun. 10 2002 


Abell 1795 (3) 




5289 


ACIS-I 


VFAINT 


9.6 


Jan. 14 2004 


Abell 2029 (1) 


0.0779 


891 


ACIS-S 


FAINT 


18.7 


Apr. 12 2000 


Abell 2029 (2) 




4977 


ACIS-S 


FAINT 


74.8 


Jan. 8 2004 


Abell 2029 (3) 




6101 


ACIS-I 


VFAINT 


8.2 


Dec. 17 2004 


Abell 478 (1) 


0.088 


1669 


ACIS-S 


FAINT 


39.9 


Jan. 27 2001 


Abell 478 (2) 




6102 


ACIS-I 


VFAINT 


9.2 


Sep. 13 2004 


Abell 2204 (1) 


0.152 


499 


ACIS-S 


FAINT 


10.0 


Jul. 29 2000 


Abell 2204 (2) 




6104 


ACIS-I 


VFAINT 


9.6 


Sep. 20 2004 



includes Abell 665, 2163, 2219, 2255, 2319, 2744, and the 
Bullet Cluster IE 0657-56. 

As a control, we have also observed five comparably 
X-ray luminous clusters that are known not to contain ra- 
dio halos: AbeU 478, 576, 1795, 2029 and 2204 (Feretti & 
Giovannini 2007). Most are highly relaxed and contain cen- 
tral cooling cores, although Abell 576 exhibits a more active 
dynamical state {eg. Dupke et al. 2007; Kempner & David 
2004a). 



uctsU Bad pixels were removed and standard grade selec- 
tions applied. Where possible, the extra information avail- 
able in VFAINT mode was used to improve the rejection 
of cosmic ray events. The data were cleaned to remove pe- 
riods of anomalously high background using the standard 
energy ranges and binning methods recommended by the 
Chandra X-ray Center. The net exposure times after clean- 
ing are summarized in Table [T] Separate photon-weighted 
response matrices and effective area files were constructed 
for each region analyzed. 



2.2 Data Reduction 

The Chandra observations were carried out using the Ad- 
vanced CCD Imaging Spectrometer (ACTS) between March 
2000 and June 2007. The standard level- 1 event lists pro- 
duced by the Chandra pipeline processing were reprocessed 
using the CI AO (version 4.1.1) software package, including 
the appropriate gain maps and updated calibration prod- 



^ Our analysis uses the latest Chandra calibration products re- 
leased in January 2009 as incorporated into CALDB version 
4.1.1. The results are robust to residual systematic uncertainties 
in the calibration (see section 3.1). 
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2.3 Spatially-resolved spectroscopy 

2.3.1 Contour binning 

The individual regions for the spectral analysis were deter- 
mined using the contour binning method of Sanders (2006), 
which groups neighboring pixels of similar surface bright- 
ness until a desired signal-to-noise threshold is met. For 
data of the quality discussed here, regions are frequently 
small enough that the X-ray emission from each can be ap- 
proximated usefully by a single temperature plasma model. 
We consider only those regions of the clusters within which 
the spectra are determined reliably, with conservative al- 
lowances for systematic uncertainties in the background sub- 
traction. In detail, we set the outermost radius to be the ra- 
dius at which the ratio of background-to-total counts in the 
0.8—7.0 keV energy range rises to 20 per cent. At this surface 
brightness threshold, a 10 per cent shift in the background 
normalization will result in an approximately 10 per cent 
shift in the temperature of an isothermal plasma model, for 
both the ACIS-S and ACIS-I detectors. The total number of 
spectral regions per cluster ranges from three to almost 100. 
Regions vary in size from a few hundred to approximately 
ten thousand square arcsec. 

For our standard analysis described in Section 3, each 
independent spectral region contains approximately 10'' 
counts. This is sufficient to determine both the metallicity of 
the cluster gas and the photon index of a power-law compo- 
nent, where present, to approximately 10 per cent accuracy. 
The temperature and density of the cluster gas are always 
determined to higher precision. 

In order to generate the detailed thermodynamic maps 
of temperature, pressure and entropy discussed in Section 
4, we have also carried out a second analysis using a finer 
binning of ~ 3,000 counts/region. For this second analysis, 
no variable power-law emission components were included 
and the metallicity was fixed at the appropriate emission- 
weighted mean values determined from the ~ lO' counts- 
per-bin results. The analysis with the finer regions is suffi- 
cient to determine the temperature to ~10 per cent accuracy 
for 10 keV gas, or ~ 3 per cent accuracy for the coolest re- 
gions in the cores of the dynamically relaxed, non-radio halo 
clusters. 

2.3.2 Background modelling 

Searching for non-thermal emission signatures in Chandra 
spectra requires careful modelling of the astrophysical and 
instrumental background emission. We have taken care to 
account for spatial and temporal variations in these compo- 
nents. 

Background spectra for the appropriate detector regions 
were extracted from the blank-sky fields available from the 
Chandra X-ray Center. These were normalized by the ratio 
of the observed and blank-sky count rates in the 9.5-12 keV 
band (the statistical uncertainty in the observed 9.5 — 12 
keV flux is less than 5 per cent in all cases). 

To enable further refinement of the background mod- 
els, for each ACIS-S observation, additional spectra were 
extracted from source-free regions of Chip 5; for ACIS-I ob- 
servations, additional spectra were extracted from regions of 
Chips 0-3, as far as possible from the cluster. These spectra 
were compared with the blank-sky fields in order to check 



for consistency and, in particular, to identify the presence 
of additional soft X-ray emission components along the line 
of sight, and/or scattered solar X-rays. Where excess soft 
background emission was detected, this was modelled and 
appropriate additional components included in the spectral 
analysis described below. Additional soft background emis- 
sion components were required in the analysis of the July 
2002 observation of IE 0657-56 and for each observation of 
Abell 2029, 2163, and 2204. AbeU 2029 in particular lies 
behind the North Polar Spur and has clear excess soft back- 
ground emission visible in the Chandra spectra. 



2.3.3 Basic Spectral Analysis 

The spectra have been analyzed using XSPEC (version 11.3; 
Arnaud 1996), the MEKAL plasma emission code (Kaastra 
& Mewe 1993), and the photoelectric absorption models of 
Balucinska-Church & McCammon (1992). All spectral fits 
were carried out in the 0.6 — 7.0 keV energy band. The 
extended C-statistic available in XSPEC, which allows for 
background subtraction, was used for all fitting. 

The default spectral model (model A) applied to each 
spatial region consists of a MEKAL plasma model at the 
redshift of the cluster with the absorbing column density 
along the line of sight fixed at the Galactic value determined 
from HI studies (Kalberla et al. 2005) The normalization 
and temperature in each spatial region are included as free 
parameters in the fits. For the initial analysis with 10* counts 
per bin, the metallicity is also allowed to vary independently 
in each region. 

A series of more sophisticated spectral models were also 
examined and the results compared to those obtained with 
model A: model B is the same as model A but with the 
line-of-sight absorption allowed to vary as a single, addi- 
tional free parameter (linked to a common value across all 
regions). Model C has Galactic absorption but introduces an 
additional power-law component into each spectral region: 
this approximates the spectrum of a range of possible non- 
thermal or quasi-thermal emission components. The normal- 
ization of the power-law component in each region is a sep- 
arate free parameter. The photon index of the power-law 
component is initially fixed to a canonical value of F = 2.0. 
(We later examined models in which the photon index of the 
power-law component was allowed to vary, either in unison 
across all regions or in each spatial region independently, see 
section 3.2). 

Finally, model D builds on model A, but introduces a 
second MEKAL component in each region, with the tem- 
perature of the second emission component linked to vary 
in unison across all regions. For computational reasons, an 
upper limit on the temperature of this component of 80 keV 



A comparison of the line-of-sight Galactic column densities 
to the target clusters determined by the HI study of Kalberla 
et al. (2005) with the earlier work of Dickey & Lockman (1990) 
shows good overall agreement, with one notable exception: the 
Galactic column density for IE 0657-56 is 6.53 X 10^" atom cm"'^ 
in the Dickey & Lockman (1990) study, whereas it is 4.89 X 10^° 
atom cm"'^ in the later Kalberla et al. (2005) work. We have 
adopted the Kalberla et al. (2005) result as the more accurate 
value. 
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Table 2. Summary of results from the analysis of the 10'' count regions for the merging, radio halo clusters using spectral models A-D. 
AC and Au are the improvements in C-statistic and change in the number of degrees of freedom with respect to the default model (A, 
or in the case of Abell 2163, B) for each subsequent fit. The value P gives the approximate probability that such an improvement would 
occur randomly. (Low P values indicate significant improvements.) F is the photon index of the power-law components in model C. kT2 
is the temperature of the additional MEKAL component in model D. 



Model A B C D 



A2319 (10 regions) 

AC - 18.5 75.6 61.8 

Au - 1 10 11 

Nh 8.09 ^-^-0.3 8-09 8.09 

r - - 2.0 - 

fcTa - - - 80^8 

P - 1.8x10^5 4.8x10-12 4.9x10"' 



A2255 (4 regions) 

AC - 10.1 25.9 28.3 

Au - 14 5 

Nh 2.50 l-ltol 2.50 2.50 

r - - 2.0 

m - - - 8ot% 

P - 1.0 X 10-^ 3.6 X 10-5 3.5 X lO"'^ 



A665 (3 regions) 

AC - 21.5 25.5 14.3 

Au - 13 4 

Nh 4.33 0.6+°-^ 4.33 4.33 

r - - 2.0 

kT2 - - - 18lf 

P - 3.9x10-6 1.4x10-5 7.0x10-3 



A2163 (17 regions) 

AC - - 54.1 22.2 

Ai/ - - 17 18 

Nh - 15-0l°;i 18-3^?:o 15.5+°; 

r - - 2.0 - 

feTa - - - 60^20 

P - - 1.0 X 10-5 0.224 



was imposed. Model D also has the column density fixed at 
the Galactic value. 



Due to degeneracies between parameters for models C 
and D, we constrain the metallicity in each region to remain 
within ±50 per cent of the values determined with model 
A. By doing this and keeping the absorbing column density 
fixed at the Galactic value, we limit the impact of unphysical 
model solutions. 



For Abell 478 and Abell 2163, the line-of-sight absorb- 
ing column density is known to differ substantially from the 
radio HI values of 1.5 x lO^'atom cm" "^and 1.09 x lO^'atom 
cm~^, respectively (Johnstone et al. 1992; Allen et al. 1993; 
Pointecouteau et al. 2004; Elbaz et al. 1995). For these clus- 
ters, model B is therefore adopted as the default model and 
models C and D also have the column density included as a 
free parameter. We caution that for Abell 478, the column 
density is also known to vary as a function of position (Allen 
et al. 1993). For this reason and for this cluster alone, the 
column density was allowed to vary independently in each 
region studied. 



3 RESULTS 

3.1 The statistical evidence for non-thermal-like 
components 

Tables [2] and |3] summarize the results from the fits with spec- 
tral models A-D for the radio halo clusters and the control 
sample of non-radio halo systems, respectively. The tables 
list the total decrease in C-statistic (AC) summed across 
all regions in each cluster, and the increase in the number 
of degrees of freedom (Ai/) for each model with respect to 
the default model (typically model A, but model B for Abell 
478 and 2163), together with the approximate probability, 
P, that such an improvement could be obtained randomlyrl 
The tables also list the best-fitting absorbing column density 
for cases where this is included as a free parameter. 

Comparing the results for the simplest single- 

To estimate the probability that the improvements to the fits 
could be obtained by random chance, we have calculated a pseudo 
F— statistic, formed by the ratio F = AC/Av. This is similar 
to the usual F— statistic formed with the chi-square estimator 
{eg. Bevington 1969), and is expected to provide useful estimates 
of the probability for cases where the spectral model provides a 
reasonable description of the data. 
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Table 2. Continued. 



Model 


A 


B 


C 


D 


A2219 


(8 regions) 








AC 


_ 


0.02 


3.3 


3.9 


Ai^ 


_ 


1 


8 


9 


Nh 


1.76 


^■''-0.3 


1.76 


1.76 


r 


- 




2.0 


- 


kT2 


- 


- 


- 


4+° '* 


p 


- 


0.89 


0.91 


0.92 


IE 0657-56 


(55 regions) 








AC 


_ 


10.3 


107.6 


88.8 


/^u 


_ 


1 


55 


56 


Nh 


4.89 




4.89 


4.89 


r 


_ 




2.0 


_ 


kT2 


- 


- 






P 


- 


1.3 X 10-3 


2.9 X 10-5 


3.4 X 10-3 


A2744 


(6 regions) 








AC 




0.3 


7.1 


17.0 


/^u 




1 


6 


7 


Nh 


1.39 


1 4+0-4 


1.39 


1.39 


T 






2.0 




kT2 








22+_f 


P 




0.58 


0.31 


0.017 



E 



I 64 
-B 

N 

I 62 
I 

a SO 




5* 

L.-. .1 t-t-J .IuIjI-IJI ■I.IjL|J.II..JI.|J|.IJI H-IJI-IjI H.!.!!..^ 



2 J 4 S 6 ? 



Figure 1. The 68 and 95 per cent confidence contours of the 
power-law and thermal normalizations summed across all regions 
in Abell 2319. Contours were determined from a Markov Chain 
Monte Carlo analysis and confirm that the power-law component 
is detected (normalization > 0) at high significance. 



temperature models, A and B, we see that for four out of 
six radio halo clusters for which both models were examined, 
allowing the absorbing column to fit freely leads to a signifi- 
cant improvement in the fit (AC > 10 for 1 additional degree 
of freedom). The most significant improvement, AC ~ 22, 
is obtained for Abell 665. In all four of these six cases, how- 
ever, the best-fit column density is significantly lower than 
the Galactic value determined from HI studies, a result that 



is physically implausible and which suggests that model B 
does not provide a complete description of the spectra. (For 
column densities < 8 x 10'^'' atom cm-^, the Galactic HI 
values are expected to be accurate to better than 10 per 
cent). The only two radio halo clusters for which significant 
improvements (>99 per cent confidence) to the fits were not 
obtained using model B over model A are Abell 2219 and 
Abell 2744. For the seventh radio halo cluster, AbeU 2163, 
which lies at low Galactic latitude, an accurate Galactic col- 
umn density estimate is not available from HI studies and 
model A was therefore not examined; model B is taken as 
the default model for this cluster. 

Importantly, none of the control sample of non-radio 
halo clusters show a comparably significant improvement 
when using model B over model Arl As discussed above, 
Abell 478 is excluded from this test since it lies at low Galac- 
tic latitude and the estimate of the line-of-sight column den- 
sity from HI studies is poorly determined. The most statis- 
tically significant detection is for Abell 2029, although the 
background model for this cluster is complex with relatively 
high systematic uncertainties (Section 2.3.2). It is unlikely 
that our modelling of the excess soft emission along the line 
of sight, associated with the North Polar Spur, is correct 
in detail (we assume it to have constant surface brightness 
across the field). Note also that the data for this cluster ex- 
hibit an exceptionally high signal to noise ratio and the best 
fit column density is only lower than the Galactic value by 
~ lO^'' atom cm-^ 

Comparing the results for models C and A, we see 
that the introduction of a power-law component with pho- 
ton index F = 2.0 leads to comparable or even larger im- 



5 We have confirmed that consistent results are obtained whether 
the ACIS-S and ACIS-I data for individual clusters are analyzed 
separately or together. 
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Table 3. Details as for table [3] but for the relaxed, non-radio halo clusters. For Abell 478, the line of sight Galactic column density is 
known to significantly exceed the value determined from HI studies (Johnstone et al. 1992; Allen et al. 1993) and vary spatially across the 
cluster. For this cluster, the Afjj value was allowed to vary independently in each region. Low P values indicate significant improvements. 
For model D, only regions beyond a radius of lOOh^^ kpc were used, to limit the impact of projection effects in the presence of strong 
temperature gradients (Section 3.2). 



Model 



D 



A576 


(4 regions) 






(4 regions) 


AC 




0.6 


7.5 


4.5 


Ai/ 


- 


1 


4 


5 


Nh 


5.50 


c c + O.S 
''■^-0.4 


5.50 


5.50 


r 


- 


- 


2.0 


- 




- 


- 


- 




P 


- 


0.44 


0.11 


0.48 


A1795 


( 24 regions) 






(8 regions) 


AC 




0.08 


31.4 


24.5 






1 


24 


9 


Nh 


1.19 


1 10+0.09 


1.19 


1.19 


r 






2.0 














P 




0.77 


0.14 


3.7 X 10-7 


A2029 


(90 regions) 






(40 regions) 


AC 




9.1 


113.0 


68.8 


Ai^ 




1 


90 


41 


Nh 


3.26 


q r,q+0.06 

'-'•'-'•^-0.07 


3.26 


3.26 


V 






2.0 




kT2 










P 




2.6 X 10-^ 


0.051 


4.2 X 10-3 



A478 


(37 regions) 






(21 regions) 


AC 






45.7 


28.3 


Au 






37 


22 


Nh 




24.7 - 31.1 


24.7- 35.8 


24.6 - 31.3 


r 






2.0 




fcTz 








41+39 


P 






0.15 


0.17 


A2204 


(7 regions) 






(3 regions) 


AC 




6.6 


1.5 


5.2 


Ai^ 




1 


7 


4 


Nh 


5.67 




5.67 


5.67 


r 






2.0 












3 1+'^-^ 

■^■^-1.2 


P 




0.01 


0.98 


0.27 



provements in the fits to the same four radio halo clusters 
{i.e. those that fitted to sub-Galactic column densities with 
model B), as well as for Abell 2163, for which a reliable 
Galactic column density measurement is not available. For 
Abell 2219, 2744, and for all of the relaxed non-radio halo 
clusters, however, no comparably significant improvement 
in the fits is obtained with model C over model A. Figure [l] 
shows the combined power-law normalization versus thermal 
normalization (in units of 10"'^) summed across all regions 
in Abell 2319. Power-law emission is detected at high sig- 
nificance. The power-law normalization is approximately 10 
per cent of the thermal normalization and, thus also, ap- 
proximately 10 per cent of the total fiux. 

We conclude, based on the improvements obtained with 
models B and C, that 5/7 radio halo clusters show evidence 
for additional model components, over and above those ac- 



counted for by the default thermal model. |j The additional 
emission components can be approximated by simple power- 
law models, although we stress that these models are un- 
likely to provide a complete, physical description of the data. 
In contrast, the relaxed, non-radio halo systems do not show 
comparable evidence for such power-law components. 

Similar conclusions on the presence or otherwise of non- 
thermal-like emission in the radio halo clusters are also 
drawn from the fits with model D, where a second MEKAL 
component, rather than a power-law model, is incorporated 
into the modelling (other than for Abell 2163; although this 
cluster also has the column density included as a free param- 

^ Similar conclusions regarding the presence or otherwise of 
power-law components in the radio halo clusters are obtained 
when the photon index of the power-law component is also in- 
cluded as a single free parameter. 
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Figure 2. (a) left panel: Spatial map of the surface brightness (in ergcm~^s~^arcsec'^) of non-thermal-like X-ray emission in the Bullet 
Cluster, IE 0657-56 {z = 0.297). A power-law component is only included in regions where it is statistically required at greater than 
90 per cent significance (Section \i.2\ . The 1.34 GHz radio surface brightness contours from Liang et al. (2000) are overlaid in black. 
(Radio point sources have been removed), (b) right panel: Map of the photon index of the power-law components, with the radio surface 
brightness contours overlaid. Only regions with 'non-thermal' surface brightness greater than 10~^® ergs s^^ cm^^ arcsec"^ are shown 
in this map. Each spatial region has ^ 10* net counts in the 0.6 — 7.0 keV Chandra band. 



eter). The nominal temperatures of the additional emission 
components are typically high, with kT2 > 15 keV. 



As discussed in Section 3.2, the majority of the con- 
trol sample of relaxed, non-radio halo clusters (all, excluding 
Abell 576) are known to contain strong temperature gradi- 
ents in their (cooling) cores. The presence of such gradients 
will naturally lead to requirements for additional tempera- 
ture components in fits to projected spectra for the central 
regions of these systems. To mitigate this effect, for the anal- 
ysis of the control sample of clusters with model D, we have 
excluded the inner lOO/ifg^kpc radius. (The reduced number 
of regions studied are indicated in column 4 of Table 3.) Of 
the control sample, Abell 2029 presents a mild requirement 
for additional temperature components, and Abell 1795 a 
more significant requirement, with most of the signal aris- 
ing from within radii ~ 200/i.fQ^kpc in both cases. Possible 
origins for these requirements are discussed in Section 5. 



Finally, to gauge the level at which residual instrument 
calibration issues may affect the results, we have repeated 
the analysis in three other energy bands: 0.6 — 4.0 keV, 
0.6 — 5.0 keV and 2.0 — 7.0 keV. For the merging, radio 
halo clusters, similar detections of additional emission com- 
ponents are made using each energy band. The control sam- 
ple of regular, non-radio halo clusters provides comparable 
null results in all energy bands studied. 



3.2 Detailed spatial mapping of the 'non-thermal' 
components 

The results of Section 3.1 provide a simple, statistical mea- 
sure of the overall significance of non-thermal-like emission 
signatures in the clusters. However, even where these overall 
signatures are strong, as in 5/7 of the radio halo clusters, we 
do not necessarily expect power-law emission components to 
be present in every region in the clusters. We have therefore 
attempted to map the spatial variation of this emission. 

Starting with model A, a single temperature model with 
the absorption fixed at the Galactic HI value, we have de- 
termined, for each individual region of each cluster, the sta- 
tistical significance of the improvement to the fit obtained 
by introducing a power-law component with a photon index 
that is allowed to fit freely in that region. Where the signif- 
icance of the improvement was found to exceed a threshold 
value (we adopt 90 per cent) the power-law component was 
kept in the model for that region]^ Where it was not, model 
A (single, thermal emission component only) was used as the 
appropriate region model. In this way, the number of regions 
requiring power-law emission components in each cluster 
were determined. Finally, the overall absorbing column den- 
sity was allowed to vary and checked for consistency against 
Galactic HI survey results]^ Maps of the surface brightness 

^ We have constrained the metallicity to lie within ±50 per cent 
of the value determined with model A. Similar results are also 
obtained if the metallicity is simply constrained to be less than 
solar. 

* An alternative algorithm to identify regions requiring power- 
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Table 4. Results from the spatial mapping of the power-law components in the radio halo clusters (Section |3.2[ ). Column 2 lists the 
Galactic column density for each cluster determined from HI studies (Kalberla et al. 2005; for Abell 2163 a reliable column density 
measurement is not available from radio data). Columns 3 lists the total number of regions for which a power-law component is detected 
at 90 per cent significance or greater. Columns 4 — 7 list the column density (in units of lO^^atom cm~^) determined from the fits, the 
mean photon index of the required power-law components, and the summed 0.6 — 7.0 keV flux (in units of 10~^^ ergs cm~^) of these 
components, respectively. (Regions with power-law photon indices, F < 1, typically exhibit very low non-thermal surface brightness. We 
consider such results to be unphysical and have excluded them from the average values in columns 5,6). 



Cluster A''H,Gal # regions A^h rmoan PL Flux 

Abell 2319 8.09 3/10 8.3^0.4 2.0I0 3 ^^^s 

Abell 2255 2.50 2/4 3.01^0° l-68to 2? ^■'^-W 

AbeU 665 4.33 2/3 2.6ti;5 l-63to.2i 4.2l^;2 

Abell 2163 - 4/17 15.4t;^;* l-Slto'os 3.9i;t;° 

Abell 2219 1.76 0/8 l-7^o| 

Abell 2744 1.39 3/6 ^■'itoA l-66to.i3 O-^laa 

IE 0657-56 4.89 13/55 4.99j:o;l4 l-SOtopj 0.95to:t5 



of the power-law emission, where detected, and the photon 
index can then be generated in a straightforward manner. 
By including 'non-thermal' emission components only where 
they are statistically required, we minimize the impact of 
random errors and residual calibration uncertainties on the 
maps. 

The best fit absorption column densities and mean 
emission-weighted photon indices determined from the de- 
tailed mapping at the 90 per cent threshold are summarized 
m Table |4] In all cases, the measured column densities are 
consistent with the Galactic values (Kalberla et al. 2005) 
at the ~ 95 per cent (2cr) confidence level. This contrasts 
with the results obtained with model B (single temperature 
model) where column densities less than the Galactic value 
was observed in most cases (Section 3.1)Q 

Excluding regions with unphysically low photon indices 
I.e. r < 1, which are almost always detected at very low 
surface brightness and which can be plausibly associated 
with residual calibration and background subtraction errors, 
the mean emission-weighted photon indices for the clusters 
range from 1.5 < T < 2.0. 

The results for the Bullet Cluster, 1E0657-56, which has 
the longest Chandra exposure components are shown in Fig 
[2^. Also shown overlaid in Fig [2^ are the contours of con- 
stant radio halo brightness from the 1.34 GHz radio data of 
Liang et al. (2000; radio point sources have been removed). 
The correspondence between the radio contours and power- 
law X-ray emission map is interesting. The peaks of the radio 
halo emission and power-law X-ray emission both occur in 
the main cluster center. There is also an extension of both 
emission components, albeit at lower surface brightness, to- 



law components was also examined, wherein the column density 
was allowed to fit freely in each region, although constrained not 
to vary by more than 50 per cent from the Galactic value. This 
algorithm led to similar results on the regions requiring power-law 
emission components. 

^ We have also mapped the spatial variation of the non-thermal- 
like emission using a more restrictive threshold set at 99 per cent 
confidence. The best-fit column densities and emission-weighted 
photon indices obtained from this analysis are broadly consistent 
with those determined using the 90 per cent threshold, although 
the column density for IE 0657-56 is reduced slightly below the 
Galactic value. 




5E-17 1E-16 1.5E-I6 2E-16 2.5E-16 3E-16 35E-1B 



Figure 3. Spatial map of the surface brightness (in units of erg 
s~^ cm~^ arcsec"^ of non-thermal-like X-ray emission in Abell 
2029. For comparison with Fig. 1, the scale has been adjusted to 
represent what would been seen if this cluster were at the same 
redshift as IE 0657-56 and the observation time were extended to 
maintain the current signal-to-noise ratio. The inner circle marks 
a radius of lOOh^^ kpc, within which the observed temperature 
declines steeply and projection effects are strong. The outer circle 
shows marks the outer limit of the analysis. 

ward the merging 'bullet' subcluster. We conclude that there 
is a correspondence between the locations of the strongest 
power-law X-ray emission and strongest non-thermal radio 
emission in the cluster. 

Fig[2|3 shows the spatial variation of the photon index of 
the power-law X-ray emission in IE 0657-56, again with the 
radio contours overlaid. The regions of strongest power-law 
fiux have a photon index F ~ 1.4 — 1.8. Regions with pho- 
ton index F < 1.0, all of which are detected at relatively low 
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surface brightness, SBpL < 1 x 10~ , are probably unphys- 
ical. Summing the total flux of all power-law components 
detected at a confidence level of 90 per cent or greater and 
with a photon index F > 1, we obtain a total 0.6 — 7.0 keV 
power law flux from IE 0657-56 of 0.95to;iJ x 10"^^ ergs s"^ 
cm~^ and an emission-weighted photon index of 1.501;q q5. 
Extrapolating this result to the 20 — 100 keV band, we pre- 
dict a flux of 3.8^0 9 X 10~^^ ergs s~^ cm"'^, consistent with 
the result of Petrosian et al. (2006) of 5 ± 3 x 10"^^ ergs s"^ 
cm"^ based on RXTE and XMM-Newton X-ray data. 

We have applied a similar mapping analysis to our con- 
trol sample of non-radio halo clusters. Fig [3] shows the re- 
sults for Abell 2029, which has a deep exposure and the most 
accumulated X-ray counts of any cluster in the study. We 
have scaled the surface brightness of the power-law emission 
in Fig. |3] to the levels that would be detected if the cluster 
were observed at the same redshift as IE 0657-56 (assum- 
ing an ultra-deep observation with sufficient exposure time 
to maintain the current signal-to-noise ratio). For relaxed 
clusters like Abell 2029 at relatively low redshifts, it is im- 
portant to recall that the temperature profiles are commonly 
observed to decline steeply in their cores (e.g. Allen, Schmidt 
& Fabian 2001; Vikhlinin et al. 2005). The analysis of pro- 
jected spectra in the innermost regions of such clusters will 
therefore always, given good enough data, provide signatures 
of multiphase gas, simply due to projection effects. The inner, 
green circle in Fig. [s] marks a radius of 100/ifg^kpc, within 
which the temperature profile in Abell 2029 is observed to 
decline steeply. The regions with the most significant detec- 
tions of power-law components lie within or near this radius; 
the photon indices of these components are typically ffat, 
with r ~ 1.0 - 1.5. (See also Section 5.2). 

The results for the other, relaxed non-radio halo clusters 
show no convincing evidence for power-law X-ray emission 
components beyond their inner regions. 



4 THERMODYNAMIC MAPPING OF 
IE 0657-56 AND ABELL 2029 

IE 0657-56 and Abell 2029 have the deepest Chandra expo- 
sures of the radio halo clusters and control sample of non- 
radio halo clusters, respectively, and are therefore the sys- 
tems best-suited to detailed thermodynamic mapping. The 
Chandra X-ray images and maps of temperature, pressure 
and entropy for the two clusters are shown in Figs. |4] and [5] 
The thermodynamic maps for IE 0657-56 are complex 
and confirm that the cluster is far from hydrostatic equi- 
librium, having experienced a recent, massive merger event. 
The peak of the thermal X-ray surface brightness is cen- 
tered on the fast-moving, merging subcluster i.e. the 'bullet' 
(see Markevitch et al. 2002; Markevitch 2006; Markevitch 
& Vikhlinin 2007 and references therein), rather than the 
overall X-ray center. The merging subcluster also contains 
the lowest temperature and lowest entropy gas. However, 
the regions of highest thermal pressure in the cluster are 
not coincident with the X-ray surface brightness peak, but 
are rather found near the overall X-ray center. The multi- 
ple apparent pressure peaks may be sites of recent and/or 
ongoing shock activity and might be expected to be loca- 
tions of ongoing particle acceleration. Fig. [Ojshows the 1.34 
GHz radio contours overlaid on the same pressure map. The 



figure confirms that the sites of highest thermal pressure 
and strongest radio emission are closely correlated (see also 
Govoni et al. 2004). 

The shock front ahead of the merging subcluster in 
IE 0657-56 has been studied by Markevitch et al. (2002), 
Markevitch (2006), and Markevitch & Vikhhnin (2007). This 
shock is also clearly visible in the thermodynamic maps of 
Fig. [4] Our analysis of the temperature jump at the shock 
front indicates a Mach number M = 2.51q | which corre- 
sponds to a velocity of « = 4000lgQQ km s~^, consistent 
with previous work. We note a hint of a second, possible 
front leading the cool subcluster in Fig. [IJj, and trailing the 
main shock (see also Markevitch & Vikhlinin 2007). 

In contrast, for Abell 2029, the thermodynamic maps, 
and especially the pressure map, appear remarkably sym- 
metric and regular. The cluster is not in per feet equilibrium, 
as evidenced by detailed structure observed in the X-ray 
image, especially at small radii (see eg. Clarke, Blanton & 
Sarazin 2004), as well as the temperature and entropy maps 
shown here (an extension of low entropy gas is observed 3 — 4 
arcmin to the southwest). However, the deviations from el- 
liptical symmetry in the pressure map are small and argue 
that the hydrostatic assumption employed in X-ray mass 
analyses is likely to be a reasonable approximation in the 
observed region of the cluster, confirming earlier conclusions 
drawn on the basis of the X-ray surface brightness distribu- 
tion alone {eg. Buote & Tsai 1996; Allen et al. 2008). The 
peak of the X-ray surface brightness coincides precisely with 
the location of the dominant cluster galaxy, IC 1101. The X- 
ray peak also marks the location of the lowest temperature 
and lowest entropy gas in the cluster, and of the highest 
thermal pressure, a configuration that is convectively sta- 
ble. 

Qualitatively similar results on the thermodynamic 
properties are obtained if power-law components consistent 
with those discussed in Section 3.2 are accounted for explic- 
itly in the analysis (although the measured temperatures 
are, typically, slightly lower). Results from thermodynamic 
mapping of other clusters in the sample are presented in the 
Appendix. 

5 THE ORIGIN OF THE POWER-LAW 
COMPONENTS 

5.1 Summary of Results 

Our analysis has revealed the presence of significant non- 
thermal-like X-ray emission in five of a sample of seven mas- 
sive, merging clusters with luminous radio halos. Our results 
shed new light on the origin of this emission. The most im- 
portant, new facts are 1) the presence of such components 
in a large fraction of the radio halo clusters studied; 2) the 
absence (or lower detection rate) of such components, at sim- 
ilar levels, in the control sample of non-radio halo clusters; 
3) the observed spectra of the non-thermal-like components, 
which are consistent with simple power-law models with rel- 
atively ffat photon indices, 1.5 < 7 < 2.0, or bremsstrahlung 
emission from hot, kT > 15 keV, gas; and 4) the spatial cor- 
relation between the regions of the brightest non-thermal 
radio halo emission, brightest, thermal X-ray emission, and 
strongest non-thermal-like X-ray signatures in the central 
regions of IE 0657-56 (Fig. 1). 
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Figure 4. Thermodynamic maps for the Bullet Cluster, IE 0657-56. Regions follow lines of constant surface brightness and have ~ 3000 
net counts, leading to la fractional uncertainties in the plotted quantities of < 10 per cent. Upper left panel: surface brightness (counts) 
in the 0.8 to 7.0 keV band. Upper right panel: temperature, kT, in units of keV determined with spectral model A. Lower left panel: 
pressure, P, in units of keV cm 2 arcsec . Lower right panel: entropy, S, in units of keV cms arcseca . The shock front analyzed by 
Markevitch et al. 2002 (see also Markevitch 2006 and Markevitch and Vikhlinin 2007) is marked by the arrow on the temperature map. 



Summarizing the statistical evidence from Section 3, we 
see that the comparison of results from models B and A pro- 
vides evidence (P < 10~^) for excess soft emission (i.e. less- 
than-Galactic absorbing column densities) in 4/6 radio halo 
clusters for which this test is possible. Of the control sam- 
ple, only Abell 2029 shows mildly significant evidence for a 
soft excess, which can be plausibly associated with residual 



uncertainties in the background modelling for this system. 
Comparing results for models C and A, we see that 5/7 radio 
halo clusters show evidence for power law components (in- 
cluding the same 4 noted above) whereas, strikingly, 0/5 of 
the control sample do. The more detailed mapping discussed 
in Section 3.2 shows that once power law components have 
been included, where statistically required, the absorbing 
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Figure 5. Thermodynamic maps for the relaxed, non-radio halo cluster Abell 2029. Other details as for Fig|4] 



column densities for the radio halo systems become consis- 
tent with the Galactic values determined from HI studies. 
The results from the tests with model D are somewhat more 
ambiguous, as might be expected given the additional free- 
dom of this spectral model: 4/7 radio halos show evidence 
for additional temperature components, as do 2/5 of the 
control sample. 



5.2 Physical Explanations 

It has been common in the literature to interpret limits on 
non-thermal X-ray emission from galaxy clusters in terms of 
simple IC models, wherein the same population of relativis- 
tic electrons responsible for the observed radio halo emission 
is assumed to Compton scatter CMB photons up into the 
X-ray band. In this specific scenario, and in the absence of 
other non-thermal or quasi-thermal X-ray emission compo- 
nents, the combination of radio and IC X-ray flux measure- 
ments can be used to estimate the magnetic field strength 
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Figure 6. Pressure map (in keV cm^2 aresec"^ for IE 0657- 
56 witli tlie 1.34 GHz radio lialo contours overlaid. Tiie strongest 
radio halo emission coincides witli the regions of highest apparent 
thermal pressure and strongest 'non-thermal' X-ray emission. 

{eg. Rephaeli 1979; Fusco-Femiano et al. 1999; Blasi & Co- 
lafrancesco 1999; Eckert et al. 2008). We have carried out 
such an analysis for the Bullet Cluster, IE 0657-56 assuming 
that the observed power-law X-ray components are entirely 
due to IC emission. In this case, 

Fx _ Eph f jx\ ^ (7Jf) _ EcMB f Jx\ ^ ^" 

Fr ^ BySn \'yR j TV (7^) " BySn {-yR J ' ^ ' 

where A*' (7) is the number density of electrons with Lorentz 
factor 7, 7x = (vxIvcmb)^^'^ , 7fl = [vR/v^ycY^'^ , vx is the 
frequency of the X-rays, vcmb is the frequency of the up- 
scattered CMB photons, vr is the frequency of the measured 
radio flux, and v^yc is the cyclotron frequency. Ecmb is the 
energy density of CMB photons at the cluster and q = F — 1 
is the spectral index of the power-law components. 

Within the context of this specific model, we infer a 
peak B field of 0.15 \iG associated with the regions of 
strongest non-thermal-like X-ray emission. Here we have 
used the photon index for each region as measured from 
the X-ray data. If we had instead assumed a fixed photon 
index F = 2.0, then the implied B fields would have been 
several times higher. 

Two factors argue against this simplest IC interpre- 
tation. Firstly, the photon indices of the power-law X-ray 
components are flatter than typical radio halo synchrotron 
spectra (2.0 ^ F < 2.3). This implies, within the context 
of the IC model, that the same population of electrons are 
not responsible for both the radio and X-ray emission. Sec- 
ondly, such low magnetic field values and spatial field con- 
figurations appear somewhat at odds with Faraday rotation 
studies that measure fields of a few, to a few tens, of /iG in 
clusters, with the strongest fields observed in cluster cores 



{eg. Kim et al. 1991; Feretti et al. 1995; Taylor et al. 1999, 
2001, 2002, 2006, Clarke, Kronberg & Bohringer 2001; Allen 
et al. 2001; Govoni & Feretti 2004). 

Theoretically, the spectrum of electrons with Lorentz 
factor 7 ~ 1000 responsible for IC X-ray emission can be 
fiatter than those responsible for the GHz radio emission. 
Indeed, this appears to be a general prediction of more de- 
tailed calculations of the coupled, time-dependent properties 
of electrons, hadrons and plasma waves in a turbulent intr- 
acluster medium, such as will be generated following a ma- 
jor merger event {eg. Brunetti et al. 2004; Brunetti & Blasi 
2005; Petrosian & Bykov 2008). A number of authors have 
also argued that the discrepancies between Faraday rota- 
tion and IC estimates of the magnetic field may be recon- 
ciled (or at least reduced) by employing more realistic dis- 
tributions for the non-thermal particles and magnetic fields 
{eg. Petrosian 2001; Carilh & Taylor 2002; Newman, New- 
man & Rephaeh 2002; Kuo et at 2003; Brunetti et al. 2004; 
Govoni & Feretti 2004). Thus, although the simplest IC sce- 
nario described by equation 1 is challenged by our results, 
more sophisticated treatments may prove consistent with 
the data. 

A second, alternative interpretation for the observed 
power-law X-ray components is bremsstrahlung emission 
from supra-thermal electrons, energized by collisions with 
an energetic non-thermal proton population in the clusters 
(Liang et al. 2002; Dogiel et al. 2007; Wolfe & Meha 2008; 
see also Ensslin et al. 1999, Sarazin & Kempner 2000, Blasi 
2000). The observed spectra appear consistent with this 
model, although physical challenges remain {eg. Petrosian 
& Bykov 2008; Petrosian & East 2008). 

Thirdly, it remains possible that a contribution to 
the observed non-thermal-like X-ray flux could arise from 
synchrotron emission from ultra-relativistic electrons and 
positrons, which would also be characterized by hard power- 
law spectra (F ~ 1.5; Inoue et al. 2005). In this case, 
however, one might expect the hard X-ray emission to be 
brighter in the outskirts of the clusters, where the accretion 
shocks are strongest. 

5.3 The potential impact of temperature and 
metallicity substructure. 

A fourth, possible explanation for the results on non- 
thermal-like, power-law signatures is the presence of signif- 
icant temperature and/or metallicity structure in the clus- 
ters. 

In order to explore the potential impact of such effects, 
we have generated and analyzed simulated X-ray data sets 
with characteristics similar to the Chandra observations of 
Abell 2319 (similar instrumental response, exposure time, 
flux, redshift, mean temperature and absorbing Galactic col- 
umn density). As with the real data, the simulated data sets 
each have ten spectral regions with 10* counts. We incorpo- 
rate a broad range of temperature and metallicity distribu- 
tions in and across the 10 regions. In each case, the mean 
emission-weighted temperatures in the ten spatial regions 
were constrained to match the observed values for Abell 
2319. The simulated data sets were fitted with the same 
basic models as the real data (models A-D) and the im- 
provements in the fits, i.e. the AC/ A// values, obtained for 
models B,C and D with respect to model A were noted. 
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We first examined wiietiier tlie observed power-law sig- 
natures could be accounted for by hot, shocked gas with 
kT ~ 20keV. The flux in the 20keV component was allowed 
to account for up to 50 per cent of the total in each region. 
In all cases, the AC/A;/ values obtained with this model 
were much smaller than the observed values. We conclude 
that the observed power-law components are unlikely to be 
due simply to the presence of unresolved hot, shocked gas 
in the clusters. 

We next examined the potential impact of broad, tem- 
perature distributions, such as might occur if the gas in each 
spatial region were highly multiphase. A range of uni-modal 
temperature models, with normal, log-normal, skew-normal 
or T-distributions of varying widths (widths ranged from 
1-5 keV for normal, skew-normal, and T-distributions or 
0.1-0.6 for lognormal distributions) were explored. Although 
for individual regions broad temperature distributions can 
sometimes generate power-law-like signatures, in no case did 
the combined results for all 10 regions provide comparable 
AC/Az^ values for models B-D to those observed. 

All of the radio halo clusters appear to have experi- 
enced recent, massive merger events. Such events can can 
introduce large amounts of relatively low temperature gas 
into cluster cores (e.g. the 'bullet' in IE 0657-56). We have 
therefore examined the potential impact of broad, bimodal 
temperature distributions in the clusters. The gas in each 
region was assumed to contain two dominant components, 
with mean temperatures and luminosities scaled according 
to standard virial relations (Mantz et al. 2009, in prep; mass 
ratios in the range 2:1 to 4:1 were examined). Each compo- 
nent was allowed to have its own broad range of temperature 
components, again represented by normal, log-normal, skew- 
normal or T-distributions of varying widths. Although this 
model provides larger AC/ A;/ values than the above, again, 
in no case did the combined results for all 10 regions provide 
comparable AC/ A:/ improvements to the real data. 

We next examined the impact of metallicity variations 
in the clusters. Galaxy clusters are commonly observed to 
contain strong metallicity gradients, with the strongest gra- 
dients observed in relaxed, cool core clusters, although gra- 
dients also observed in some merging systems [eg. Allen 
& Fabian 1998, De Grandi & Molendi 2001, Leccardi & 
Molendi 2008). When observed in projection, inner high- 
metallicity gas will be grouped with outer, lower-metallicity 
material, which can complicate the interpretation of spectra. 
We have examined the potential impact of such effects. We 
assume that each spectral region contains a bimodal metal- 
licity distribution, with metallicities differing by factors from 
2:1 to 10:1, and relative fluxes in these components ranging 
from 1:10 to 10:1. In the flrst case, with isothermal gas 
in each region, the impact of the metallicity variations is 
limited. However, if both metallicity and temperature vari- 
ations are allowed - be it in a uni-modal or bimodal man- 
ner for the temperature - then significant AC/Ajy values, 
comparable to the observations, can be obtained. The max- 
imum AC/Ai/ values are observed for roughly equal fluxes 
in the low and high metallicity components. For such cases, 
we have also confirmed that allowing the power-law index 
in each region to vary independently leads to consistent re- 
sults on the distributions of photon indices, as well as fiux in 
the inferred power-law components. For example, a model 
with a 3:1 mass merger, a T-distribution of temperature, 



and a 1:1 mix of high and low metallicity gas satisfies all 
requirements. 

It is important to consider why, if temperature and 
metallicity variations can together explain the observed 
power-law signatures in the clusters, are these only detected 
in radio halos? A possible resolution lies in the differences 
in central X-ray surface brightness of relaxed cool core and 
merging, radio halo systems: the central surface brightness is 
typically an order of magnitude higher in relaxed cool-core 
clusters (Figure [7|. For example, cool core clusters might 
have beta-model core radii of 0.02 virial radii or less, whereas 
for massive, merging systems a value more like 0.2 virial 
radii is observed (e.g. Allen 1998). For illustration, we con- 
sider a step-like metallicity distribution ranging from 0.5 so- 
lar within 0.2 virial radii to 0.1 solar beyond. When viewed 
in projection, a region spanning the central lOOkpc radius 
for the cool core cluster will have only ~25 per cent of the 
flux in the lower metallicity component. In contrast, for the 
merging system with the larger core, the flux in the outer, 
lower metallicity component will be closer to 50 per cent, 
for which maximum AC/Az^ values with models B-D, with 
respect to model A, are observed. We note that this may 
also explain the requirements for extra spectral components 
at radii of 100 - 200/ifo^ kpc using model D in AbeU 1795 
and AbeU 2029, since for these annuli also the contribution 
of the projected to the total measured flux approaches 50 
per cent. (For AbeU 2029, however, it also remains possible 
that the signature may in part be generated by processes 
associated with the mini radio halo in this cluster reported 
by Govoni et al. 2009, which also has an extent of approxi- 
mately 200ft^g^kpc.) 

Finally, for completeness, we have also carried out sim- 
ulations in which true, non-thermal power-law components 
were included in the clusters, with photon indices in the 
range 1.25 < F < 2.5 and approximately 10 per cent of the 
total flux in the power-law component (as observed). For 
photon indices in the range 1.5 — 2.0, AC/Aj^ values consis- 
tent with the observed values are also be obtained. 

We conclude that temperature and metallicity varia- 
tions, similar to those described above, could be responsible 
for the observed, non-thermal-like, power-law components 
in the clusters. However, these signatures may also be due 
to true non-thermal or quasi-thermal processes, as discussed 
in Section 5.2. 



6 SUMMARY 

Using high-quality Chandra data, we have presented new ev- 
idence for spatially-extended, non-thermal-like X-ray emis- 
sion components in five of a sample of seven luminous, merg- 
ing, radio- halo clusters. A control sample of five, comparably 
X-ray luminous and mostly relaxed clusters without radio 
halos do not, in general, exhibit similar emission signatures. 

The spectra of the detected non-thermal-like emission 
components can be approximated by simple power-law mod- 
els with photon indices 1.5 < F < 2.0. The data for the 
Bullet Cluster IE 0657-56, demonstrate a spatial correlation 
between the regions of highest thermal pressure, strongest 
radio halo emission, and strongest power-law emission in the 
cluster core. Our detailed X-ray maps reveal complex ther- 
modynamic structure in the radio halo clusters, driven by 
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Figure 7. 3-dimensional representation of the projected surface brightness for the cool core cluster Abell 2029 (left) and the radio halo 
cluster Abell 2319 (right) scaled to appear as they would if observed at the same redshift. The flat surface brightness core of Abell 2319 
with respect to that of Abell 2029 (core radius of 120 kpc versus 20 kpc, respectively) is the most obvious morphological distinction and 
impacts on the relative importance of projection effects in the two systems. The X and Y axes span 1 Mpc on a side. The Z axis shows 
the surface brightness in units of counts arcmin"^. 



the ongoing merger activity. We confirm the presence of the 
previously-discovered shock front in IE 0657-56. We also re- 
port the discovery of a clear, large-scale shock front in Abell 
2219, as well as indications of further shock activity in other 
clusters (see Appendix). 

We have highlighted several possible explanations for 
the observed, power-law X-ray emission components. One 
possibility is IC emission; however, the spectra of the elec- 
tron populations responsible the IC X-ray emission must 
be flatter than those producing the radio halo flux. The 
magnetic field strengths {B ~ 0.15/iG) implied by the IC 
model in the regions of the strongest non-thermal-like X- 
ray emission are also less than those typically inferred from 
Faraday rotation measurements. Both difficulties can, how- 
ever, be understood (at least partially) within the context 
of more sophisticated treatments of particle acceleration in 
a turbulent intracluster medium (eg. Petrosian 2001; Carilli 
& Taylor 2002; Newman, Newman & Rephaeli 2002; Kuo 
et al. 2003; Brunetti et al. 2004; Govoni & Feretti 2004; Pet- 
rosian & Bykov 2008). A second interpretation in terms of 
bremsstrahlung emission from supra-thermal electrons, en- 
ergized by collisions with a non-thermal proton population, 
(Liang et al. 2002; Dogiel et al. 2007; Wolfe & Mefia 2008; 
see also Pfrommer 2008) also appears consistent with the 
observed X-ray spectra, although physical challenges with 
this model remain {eg. Petrosian & Bykov 2008; Petrosian 
& East 2008). A contribution to the non-thermal X-ray flux 
due to synchrotron emission from ultra-relativistic electrons 
and positrons (Inoue et al. 2005) is possible, although we 
might expect such a component to be brightest in the clus- 
ter outskirts, where the accretion shocks are strongest. 

Fourthly, and perhaps most interestingly, we have 



shown that power-law signatures, similar to those reported 
here, can also be generated by temperature and/or metallic- 
ity substructure in the cluster, particularly in the presence 
of metallicity gradients. In this case, an important distin- 
guishing characteristic between radio halo and non radio 
halo, cooling-core clusters is the more sharply-peaked X-ray 
emission of the latter population. 

Our results have clear implications for the interpre- 
tation of 'soft-excess' X-ray emission from galaxy clusters 
{eg. Durret et al. 2008 and references therein) and show that 
apparent soft excesses (or best-fit absorption of less than the 
Galactic value) can arise from inadequate modelling of the 
emission spectra in the presence of significant spatial vari- 
ations in temperature and metallicity, as well as additional 
emission components. 

Further X-ray and radio mapping, coupled with new 
7-ray and TeV observations and detailed hydrodynamical 
simulations will be vital for improving our understanding of 
the origin of non-thermal and quasi-thermal X-ray emission 
from galaxy clusters. 

Finally, we note that analyses such as ours push at the 
limits of what is possible with Chandra. However, our re- 
striction to using only data gathered above a conservative 
background-to-total counts threshold, the null results ob- 
tained for the control sample of non-radio halo clusters, the 
insensitivity of the results to the precise energy band used 
all argue that the basic findings are likely to be robust. 
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APPENDIX: DETAILED THERMODYNAMIC 
MAPPING OF OTHER CLUSTERS IN THE 

SAMPLE 

This appendix includes thermodynamic X-ray maps for 
all other clusters in the sample (in addition to IE 0657- 
56 and Abell 2029) that have ~20 or more regions 
containing 3,000 counts. Our results are broadly consis- 
tent with, but refine and expand on, previous studies 
{eg. Markevitch & Vikhlinin 2001; Markevitch, Vikhhnin & 
Mazzotta 2001; Fabian et al. 2001; WorraU & Birkinshaw 
2003; Sun et al. 2003; Kempner & David 2004a; Kemp- 
ner & David 2004b; Boschin et al. 2004; Govoni et al. 2004; 
O'Hara, Mohr & Guerrero 2004; Sanders, Fabian & Taylor 
2005; Sakelliou & Ponman 2006; Dupke et al. 2007). 

A key result is the discovery of a large-scale shock front 
in Abell 2219 (Fig. 8b). The cluster shows a clear, arc-shaped 
region of hot, kT > 16 keV (projected) gas roughly 2 ar- 
cmin to the northwest of the cluster center. This gas also has 
high pressure and entropy and is almost certainly associated 
with a merger-driven shock front. Together with IE 0657-56 
and Abell 520, Abell 2219 is only the third X-ray luminous 
cluster in which a clear, large-scale shock front, viewed ap- 
proximately edge-on, has been discovered (see Markevitch & 
Vikhlinin 2007 for a review). For more detail see Applegate 
et al. 2009, in prep. 

We have estimated the Mach number and velocity of 
the shock in Abell 2219. The deprojected pre- and post- 
shock temperatures of ISl^ keV and 25lg^ keV, respec- 
tively, lead to an estimate for the Mach rmmbcr of 1.9l[5 g. 
A density jump of a factor ^^1.5 is also seen at the front, 
which leads to an independent and consistent estimate for 
the Mach number of I.sIq j- Given an ambient sound speed 
of ~ 1850 km/s, estimated from the pre-shock temperature 
in the northwest region, we determine a shock velocity of 
~ 2500 km/s. Deeper Chandra observations are required to 
enable a more robust investigation of the shock properties 
in Abell 2219. Interestingly, we note that Abell 2219 also 
exhibits a region of exceptionally hot, kT = 22'tl kcV, high 
pressure gas coincident with the peak of the X-ray surface 
brightness, which is also, presumably, associated with recent 
shock activity in the cluster core. 

Abell 2744, like Abell 2219, docs not provide signifi- 
cant evidence for non-thermal-like X-ray emission. Its X-ray 
maps show a number of interesting features. The cluster is 
in the process of a massive, merger event. The overall pres- 
sure maocimum is coincident with the X-ray surface bright- 
ness peak and is associated with the main (southeastern) 
cluster core. The hottest gas {kT ~ 11 keV) is also found 
in this region. The complex entropy structure in the main 
cluster suggests possible shocks and significant ongoing tur- 
bulence possibly indicating that the colliding cores have al- 
ready passed through each other. The coolest {kT ^ 6 kcV), 
lowest entropy gas in Abell 2744 resides in the core of the 
merging subcluster, to the northwest of the main cluster. 

For Abell 2319 and 2163, the thermodynamic maps arc 
complex and confirm that the clusters have experienced mas- 
sive, recent merger events. For both systems, the tempera- 
ture, pressure and entropy maps show a high degree of asym- 
metry. As with IE 0657-56, we observe a close correlation 
between regions of strongest apparent dynamical activity 
and brightest radio halo emission in the clusters (see Govoni 
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et al. 2004 for the radio images). In the case of Abell 2319, 
wc observe a clear extension of cool, low entropy gas to the 
northwest; the radio halo emission is also extended in this 
direction. For Abell 2163, the coolest, lowest-entropy gas is 
located to the west of the main X-ray brightness peak. As 
with the Bullet Cluster, this gas may represent the surviv- 
ing core of a recently merged subcluster. Previous claims of 
hard, non-thermal X-ray emission from Abell 2319 and 2163, 
based on RXTE data, were reported by Gruber & Rephaeli 
(2002) and Rephaeli, Gruber & Arieli (2006), respectively. 
(An upper limit for Abell 2163 was also reported by Feretti 
et al. (2001), using BcppoSax.) These results, which probe 
harder X-ray energies but with limited spatial resolution, 
are broadly consistent with ours. 

For Abell 478 and 1795, the thermodynamic maps indi- 
cate relatively relaxed thermodynamic states with a high de- 
gree of symmetry. The peaks of the thermal X-ray brightness 
coincide with the locations of the dominant cluster galaxies. 
The X-ray peaks mark the locations of the lowest tempera- 
ture and entropy, and highest thermal pressure. Unsurpris- 
ingly, and as with Abell 2029, neither cluster is in perfect 
equilibrium, as evidenced by detailed structure observed in 
the X-ray image, especially at small radii (see eg. Fabian 
et al. 2001; Markevitch, Vikhlinin & Mazzotta 2001; Sun 
et al. 2003). In particular, for Abell 1795, some approximate 
north-south elongation is visible in the pressure and entropy 
maps that may be related to recent dynamical activity and 
the observed peculiar motion of the dominant cluster galaxy 
(Fabian et al. 2001). 
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Figure 8. Thermodynamic maps for Abell 2219. Note the clear shock front approximately 2 arcmin to the northwest of the cluster 
center. The density and temperature jumps at the shock front indicate a Mach number of 1.3-1.9 (see also Applegate et al. 2009, in 
prep). Note also the very hot gas (fcT = 22^^ keV) coincident with the X-ray surface brightness peak which may also be associated with 
recent/ongoing shock activity. Other details as for Fig.H 



© 0000 RAS, MNRAS 000, 000-000 



20 E. T. Million and S. W. Allen 




0.001 0.002 0.003 0.004 0.005 500 1000 1500 2000 2500 3000 

Figure 9. Thermodynamic maps for Abell 2744. The cluster is in the process of a major, ongoing merger event. The complex entropy 
structure in the main cluster core suggests significant dynamical activity. The coolest, lowest entropy gas is found in the core of the 
merging subcluster, to the northwest. Other details as for Fig.|4] 
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Figure 10. Thermodynamic maps for Abell 2319. The maps reveal the complex dynamical state of the cluster. We observe a clear 
extension of cool, low entropy gas to the northwest. Other details as for Fig.|4] 
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Figure 11. Thermodynamic maps for Abell 2163. The maps reveal the complex dynamical state of the cluster. The coolest, lowest- 
entropy gas is located to the west of the main X-ray brightness peak. As with the Bullet Cluster, this gas may represent the surviving core 
of a recently merged subcluster. The maps also indicate a possible shock front ~ 1.3 arcmin to the north of the main surface brightness 
peak. Other details as for Fig. |4] 
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Figure 12. Thermodynamic maps for Abell 478. The cluster exhibits a relaxed dynamical state. Other details as for Fig.|4] 
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Figure 13. Thermodynamic maps for Abell 1795. The cluster exhibits a generally relaxed dynamical state, although some asymmetry 
is visible in the pressure and entropy maps that may be related to recent dynamical activity and the observed peculiar motion of the 
dominant cluster galaxy (Fabian et al. 2001). Other details as for Fig.H] 
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